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Abstract 



We calculate the one loop effective potential for a 5D SU(2) gauge field theory at 
(^ I finite temperature T = 1//3 and radius R = 1/M. This calculation is performed, for the 

^J ■ first time, in the case of background fields with two constant components A"^ (directed 

towards the compact extra dimension with radius R) and A^ (directed towards the 
compact Euclidean time with radius /?). This model possesses two discrete symmetries 
known as Zm(2) and Zt{2). The corresponding phase diagram is presented in Ref. |^. 
^^' However the arguments which lead to this diagram are mainly qualitative. We present 

(~| ! a detailed analysis, from our point of view, for this phase diagram, and we support our 

arguments performing lattice simulations for a simple phenomenological model with 
two scalar fields interacting through the previously calculated potential. 



1 Introduction 

It has been noted long ago by G. 't Hooft in Ref. P that a pure 4D SU(N) gauge field theory 
at finite temperature T develops a global symmetry which is known as Z{N) symmetry. In 
Refs. 12 EI, the one loop effective potential in the presence of a constant background gauge 
field Aq was calculated. This result, which is reliable only in the weak coupling regime, 
implies a violation of the Z(N) symmetry. This is interpreted as the phase transition to the 
deconfining phase of an SU(N) gauge field theory, that is expected for high temperatures. 

In recent years, there has been an interest for models with extra compact dimensions. 
The simplest way to extend the above model, is to add an extra compact dimension y with 
radius R = 1/M. As a consequence, the model develops an additional Z{N) symmetry. 
To distinguish the two Z{N) symmetries of the model, we will call Zt{N) the one that 
corresponds to the compact Euclidean time and Zm{N) the one that corresponds to the 

*E-mail: kfarakos@central.ntua.gr 
tE-niail: paul@central.ntua.gr 



extra compact dimension (for details see the next section). Due to these symmetries the 
model possesses four distinct phases. 

A schematic phase diagram (see Fig. ^ below) has been presented, for the first time, by 
C P. Korthals Altes and M. Laine in Ref. |3]. However the arguments, in Ref. |1], which 
lead to this phase diagram are mainly qualitative. For this reason lattice simulations for a 
5D and 4D SU(2) gauge field theory at finite temperature and radius were performed in Ref. 
j^. The lattice results for d = 4 confirm ^ the phase diagram of Fig. ^ Even in the case 
of (i = 5 where the theory is not renormalizable, for fixed lattice spacing, the qualitative 
features of the above mentioned phase diagram are evident in the lattice results. 

In this paper we calculate in detail the one loop effective potential for a 5D SU(2) gauge 
field theory at finite temperature T and radius R. This calculation is performed in the case 
of background fields with two constant components A^ and A^. This result generalizes a 
previous calculation, in Ref. |lj, for one constant gauge field component Ay and A^ = 0. 

We will study whether it is possible to derive the qualitative features of the phase dia- 
gram of Fig. ^by using the perturbative result for the effective potential. Unfortunately the 
expected restoration of the Zm{N) symmetry ^ (when the system passes from region (A) to 
(B) in Fig. Q) for high temperature can not be established just from the effective potential. 
For this reason we construct a simple phenomenological model, that incorporates the fluc- 
tuations of the scalar fields, by adding to the effective potential kinetic terms. Numerical 
simulations on lattice for this model give a phase diagram that exhibits the main features of 
the expected phase diagram of Fig. [T] 

2 The Zt{N) x Zm{N) symmetry. 

The object of study, in this section, is an SU(N) gauge field theory, in a d-dimensional 
space-time at finite temperature, with one extra compact dimension. This extra dimension 
will be noted by y and it varies from to -R = 1/M, where M is the mass scale of Kaluza- 
Klein modes. In the case of finite temperature T we have another compact dimension, the 
Euclidean time r = it, which varies from to (3 = 1/T. 
The partition function of this model is: 

Z= f VA^e-'^I^yJo'^-I'^'-'^^-^^--^^' (1) 

Jb.c 

with the periodic boundary conditions 

A,{0,T,x) = A,{R,T,x) (2) 

A^{y,0,x) = A^{y,(3,x) (3) 

where x = (xi, X2, ...,Xd-2)- The field tensor is given by the equation Ff^^, = d^Ay — dyA^ — 
igd[A^,Ay] (z/, /i = l,2...d), where gd is the d-dimensional coupling constant. In addition, 
A^ = A'^T'' ( a = 1, 2...N^-1), and T" satisfies the commutation relation [T", T^] = i/'^^'^T^. 
Also note that for the components Ad-i and Ad of the gauge field we will use the notation 
Ar and Ay. 

"'^The notation d means that the model has d-2 noncompact and two compact dimensions. 
^In this work we analyze the N=2 case. 



We would like to emphasize that this action, if c? > 5 (in this work we will study the case 
oi d = 5), corresponds to a non renormalizable field theory. However this theory is viewed 
as an effective theory valid up to a finite cut-off A, and describes the low energy behavior of 
a fundamental renormalizable field theory, which may be a string theory. In this way all the 
observables of this model are rendered finite. We note that an observable like the one-loop 
effective potential, which is computed in the next section, is finite and cut-off independent. 
Also we note that the scale A is assumed to be much larger than the temperature T and the 
mass scale M (or T << A and M << A). An extensive discussion on this topic is presented 
by K. R. Dienes et al in Ref. [H]. 

The action of this model is invariant under Gauge transformations 

A' = UA,U^ + —Ud,U^ (4) 

Note that the transformed gauge fields A' should remain periodic, otherwise we would 
have violation of the boundary conditions (2) and (3) of the path integral. We see that 
these conditions are satisfied if gauge transformations are also periodic, namely U{0, r, x) = 
U{R, r, x) and U{y, 0, x) = U{y, (3, x). 

However the class of the gauge transformations that preserve boundary conditions (2) 
and (3) is wider. In this class we have also to include and the gauge transformations with 
the property: 

ziU{0,T,x) = U{R,T,x) (5) 

Z2Uiy,0,x) = U{y,P,x) (6) 

where Zi and Z2 are elements of Z{N). This means that this model possesses an additional 
global discrete symmetry Zt{N) x Zm{N), where Zt{N) corresponds to Euclidean time r 
and Zm{N) to the extra dimension y. 

Whether the symmetries Zt{N) and Zm{N) are violated or not, depends on two order 
parameters < Pr > and < Py > ^, where 

Prix, y) = —TrVe''" ^o '^-^Ay.r^x) ^^^ 

Py{x,T) = —TrVe'^''Io''^y^y(y'^''^^ (8) 

are the Polyakov loops in these directions. 

Performing a Gauge transformation, with properties (5) and (6), we see that the two 
order parameters are not invariant and transform as: < Pr >—>■ zl < Pr > and < Py >— > 
Z2 < Py >■ So depending on the values of T and M we have four possible distinct phases: 

(A) If < Pr >y^ and < Py >^ then both the symmetries Zt{N) and Zm{N) are 
violated, and then we can use a low energy effective theory which is characterized as 
3D SU(N)+adjoint matter. 

■^We remind readers, that < Pr > and < Py > are the mean values of Pr and Py in the corresponding 
functional integral. 
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Figure 1: The expected phase diagram for the 5D SU(N) at finite temperature and radius. 
This phase diagram was proposed in Ref . jl] , and it was confirmed by lattice simulations in 
Ref. ISj. The regions (A), (B), (C) and (D) in the figure are explained in the text below. 

(B) If < P^ >7^ and < Py >= then the symmetry Zt{N) is violated but the symmetry 

Zm{N) is not violated, and the theory is characterized as 4D SU(N)+adjoint matter. 

(C) If < Pr >= and < Py >y^ then the symmetry Zt{N) is not violated but the sym- 

metry Zm{N) is violated, and again the theory is characterized as 4D SU(N)+adjoint 
matter. 

(D) If < Pr >= and < Py >= then the symmetries Zt{N) and Zm{N) are not 
violated, and so there is no low energy effective theory description, then our theory is 
characterized as 5D SU(N) theory. 

As we see in Fig. ^ the M-T plane is separated into four regions every one of which 
corresponds to one of the above mentioned cases (A), (B), (C) and (D). 



3 One loop effective potential for 5D SU(2). 

In this section we will concentrate on the case of SU(2) for d = 5. We aim to compute the one 
loop effective potential in the presence of a background field with two constant components 
Ar and Ay which are directed toward the same direction in the group space. 

We split the gauge field into a classical background field 5^ and a quantum field a^ 



(/x 



1,2,3,4, 



5): 



A^ = B^ + a^ (9) 

The background field B^ is chosen to be zero in the case of noncompact dimensions x and 
constant for the compact dimensions r and y. We emphasize that in this work we study only 
background fields with zero classical energy (or F^y = 0). This happens only if we choose 
the gauge field components B^- and By toward the same direction (toward the generator T3) 



in the group space. This choice is also supported by the fact that it is a saddle point of the 
constraint effective potential as we show in the appendix. 

So the background field is chosen according to the following equations: 

„ 2nTu^ , , 

Br = T3 10) 

95 

By = -—^Ts (11) 

fl'5 

where we have introduced the dimensionless scalar fields q and u. 

We use a gauge-fixing condition of the form D'^cC = 0, which is known as background 
Feynman gauge, where 

D- = 5"^<9^ + (75£"''5j (12) 

is the covariant derivative in the adjoint representation. 

The lagrangian can be separated into three terms Lpp = \Tr {F^i^Y + Lqf + Lgt, where 
Lgf = ~|(-D^a^)^ is the gauge fixing term and Lgt = f]"-((—D'^)"-^ — D^e'^^^a^^ri'^ is the 
ghost field term. If we keep only the quadratic terms in the quantum fields we have: 



1 



Lot = -a"' 



QT o"M 



'-D^Y^5,, 



al + r^'^i-DTv" (13) 



Note that the linear term, in quantum fields, is identically zero for the case of the background 
field of Eqs. (10) and (11), as it is a solution of the equations of motion. 
The effective potential is defined by the equation: 

e-R^yy^ffM = Jv% Vt], V&, e-^y^^'"^^""'"^'''^^ (14) 

where V is the space volume. 

Integrating out the fluctuations ci^ and the ghost fields r^^, we obtain the following 
expression for the effective potential: 

K//(g, u) = ^TM{^-\)Tr ln(-D2) (15) 

In order to perform the trace in the group space we write —D^ in the following matrix form 

i -dl^ 9l{B^^ ^ 9l{B^y? 2g,Bldy + 2g,Bldr \ 

-D'=\ -2g,Bldy-2g,Bldr -dl + gl{Blf + gl{B^y? (16) 

V Q -dl) 

where we have used Eq. (12). 

Taking into account Eqs. (10) and (11) we can write the eigenvalues of the above matrix 
into the form 

Ai = -9^ - {dr - i2-nTuf - {dy - i2T:Mqf 
A2 = -(? - {dr + i2'KTuf - {dy + i2TxMqf 

M = -d'-dl-dl 



For the trace in the functional space we will use a plane wave basis u{x, r, y) ~ ^^p-x ^t2nTnT ^i2nMmy 
(m, n = 0, ±1, ±2, ..), then from Eq. (15), if we renormalize by subtracting the effective po- 
tential with no background field present (or Veff{q,u) —>■ Veff{q,u) — K//(0,0)) , we obtain 



+ 00 „ 

V,ff{q,u) = {d-2)TM Y. J 



^d-2p fp^ + {27rMf{m + qf + {27TTy{n + uf 



(27r)'^-2 \^ p2 + (27rM)2m2 + {2nTfn 

From the integral representation ln(a/6) = — /Q^°°((is/s)(e^"* — e^^'*), we obtain: 

In 



'p^ + (27rM)2(m + qf + {2-KTf{n + uf^ 



p2 + (27rM)2m2 + (27rT)2n2 / 

-°° rfs^_p2^ ((27rM)2(m+g)2 + (27rT)2(n+«)2)s _ g-({27rM)2m2 + (27rT)2„2)sj M^\ 

/O S 

If we perform first the integration over momentum we obtain: 

V^,M^u) = {d-2)TMj-^^^^^^j^ ^J{q,u,s) (18) 

where 

+ 00 

f{q,U,s) = - V [e'(('^^'^)Hm+q)'' + {'i^T)Hn+uf)s_^-{{2nMfm^ + {27TT)^n^)s^ ^^q^ 

n,m=— CO 

Using the Poisson formula 

+ 0O +00 „^oo 

^ F(n) = E [ / c^a;e2"^"^F(a;)] (20) 

n=-oo — -- "'-°° 

we obtain 



rJ=— oo ' ' 



Setting s = l/{AM^t) 

1 +00 /• + 00 

■ E ' 

r,l=—oo 

If we set t = zjir^ + l^ j (?\ and perform the integration over z, we obtain 



1 -pCXD /j_j_QQ ^^^^ 



V{rll2\ +°° 1 

V;,,(g,.) = (rf-2)M'^4^ E ,.2^/2/.2,./2 (l-^-^""^"^"") (23) 



TT^/' wf^oo (^' + ^Vp')'^/' 



where p = T jM^ and we have used the equation: /g^°° dz z'^l'^^^e^^ = T(d/2). 

From Eq. (23) the one loop effective potential for the two scalar fields can be put into 
the form: 

V^ffiq, u) = Vjffiq) + V^lfiu) + V:ff{q, u) (24) 

6 



where 



Kf/(g) = 4(rf-2)M'^^^'^/^) 



~oo 



Vl,j{u) 



^d/2 Esin'(vrrg)[- 
4(d-2)T^^£sin^(7rM[ 



TT 



d/2 



1=1 



,T{d/2) 



7C 



d/2 



1=1 r=l 



I 



+ 00 r, 

sin^(7r/M) sin^(7rrg) 

(r2 + /Vp2)d/2 



Note that V;//(g, 0) = VHAq) and 14// (0 



(25) 
(26) 

(27) 



We observe that the potentials Vjjjiu) and K//(Q') i^ Eqs. (25) and (26) are positive and 
the potential V^?j{q,u) in Eq. (27) is negative. One may think that the effective potential 
Veff{q, u) in Eq. (24) exhibits a local minimum for q = 1/2 and u = 1/2. However numerical 
calculation, for several values of T and M in all characteristic regions, shows that this is not 
the case. A typical plot of the effective potential is shown in Fig. |21 






s 




Figure 2: A-^(TM)- Ve//(g, u) as a function of q and u for M/A = 0.5 and T/A = 0.75. 

An interesting feature of the effective potential is that it is periodic in the dimensionless 
fields q and u with corresponding periods equal to one. For this reason we have plotted the 
Veff{q, u) only in the region [0, 1] x [0, 1]. Note that this periodicity of the effective potential 
is a consequence of the Zt(2) x Zm{2) symmetry of the gauge field theory. 

We could not find an analytical expression for the effective potential. However the ef- 
fective potential M~'^Vg//(g) (or T~'^V^jj{u)) can be approximated very well by a function 
of the form c g"(l — qY where the parameters c and a are determined by a nonlinear fit 



p 


c 


c/p 


a 





18.87 


- 


2.183 


0.5 


17.78 


35.57 


2.110 


1 


21.11 


21.11 


2.015 


2 


39.49 


19.74 


2.000 


3 


59.22 


19.74 


2.000 


4 


78.96 


19.74 


2.000 



Table 1: A very good approximation for the effective potential M ^VMAq) is given by a 



curve of the form c q°'{l — q)"-. The parameters c and a are determined with a nonlinear 
fit procedure for several values of p and are presented in the above table. We see that the 
parameter c is proportional to p for p » 1. The relative errors for the parameters a and c 
are of the order of a thousandth or smaller and are not presented in this table. 



procedure ^. In Tabled we present some values of c and a for several p for rf = 5. We see 
that as p increases the parameter c increases linearly with p and the parameter a tends to a 
constant value equal to 2. 

Now for the special case of p = we have: 



M-V,^(g) = ^ E ^^^^ = 18.87 g2-i83(i _ ^)2.i83 



TT^ 



(28) 



r=l 



For the case of p ^ +oo we have shown, performing accurate numerical computations, that: 



1 +00 



4 1 



1=1 



y.2 + /2/p2)5/2 



(29) 



Thus from Eqs. (29) and (25) 



M-v,y/g) 



12 +2^sin2(7rrg) 



TT 



r=l 



;P 



~ 27rV g^(l - g)^ 



(30) 



where we have used x^il- 



■x] 



6_ Y^+oo siii^(7rr3;) 
^4 Z^ 



=1 



(for this formula see for example Ref . ^HI ) • 



According to Fig. ^we expect a restoration of Zm{'2) symmetry above a temperature T^ 
{Tc > M) where the system passes from region (A) to region (B). However the perturbative 
results can not explain the restoration of the 2'm(2) symmetry for large temperatures, as 
the barrier that separates the vacua g = and g = 1 increases linearly with the temperature 
(this is also emphasized in Ref. |lj). 



4 A 3-dimensional model with two scalar fields. 



In this section we will present an analysis, from our point of view, for the 2'm(2) symmetry 
restoration and more generally for the phase diagram of Fig. ^ taking into account the 
result for the effective potential of Eq. (24). 

"'We have assumed that qe.[Q, 1]. Strictly we should write M^^Veff{q) « c [q mod 1]''(1 — [q mod 1])". 



We consider that the fields u = g^B^/2'KT and q = g^By/2'KM^ are not constant, as it 
was assumed, but they are dependant on the spatial coordinates x = {xi,X2,X3). Now we 
can construct a new action by adding to the one loop effective potential the kinetic terms 
which are obtained by substituting Eqs. (10) and (11) in the original action of the five 
dimensional gauge theory. Then we have: 



Sefflq, ^] = ¥ / d'^ (dq)' +'^ J d'x {duY + 7^ j d'x V,ff{q, u) (31) 



This simple model is viewed as a quantum field theory and the expectation value of an 
observable quantity 0{q,u) (an example of an observable quantity is the Polyakov loop in 
Eq. (38) below) is obtained by the path integral 

Note that this model is nonrenormalizable as the potential is periodic and thus includes 
powers of q and u larger than six ^. However, in this paper, our model is viewed as a low 
energy effective theory that is valid up to a finite momentum cut-off A. In this way all 
observable quantities are rendered finite. The momentum cut-off A could be identified with 
the momentum cut-off of the original gauge theory. 

Our purpose, for the introduction of this scalar model, is to incorporate fluctuations for 
the scalar fields q and u. An estimation of the intensity of fluctuations is given by the 
inverse of the coefficients in the kinetic terms of Eq.(31). The fluctuations for the fleld u 
are controlled by the parameter gl/ p and for the fleld q by g^p where p = T/M. So, for 
example, when we increase the temperature T keeping the mass scale M flxed we increase 
the fluctuations for q and suppress the fluctuations for u. 

Note that the model possesses four topologically nonequivalent vacua: (q = 0,u = 0), 
{q = 0,u = 1), {q = l,u = 0), {q = l,u = 1). These vacua are separated by potential 
barriers, as we see in Fig. |21 According to the above model, the system can jump from one 
vacuum to another only due to fluctuations of the dimensionless scalar flelds q and u. 

When our system is in region (A) (see Fig. P) q and u are frozen to one of the four vacua 
of the model. As the temperature T increases, and the mass scale M is kept flxed, the barrier 
between the vacua (g = 0, m = 0) and {q = 0,u = 1) (or (g = 1, u = 0) and {q = l,u = 1) ) 
increases rapidly (like T^ as we obtain from Eq. (26)). In addition the fluctuations for the 
fleld u are getting weaker, thus they can not help the system to jump from one vacuum to 
another and the fleld u remains frozen to m = or to u = 1. 

On the other hand, the barrier between the vacua {q = 0,u = 0) and (q = l,^ = 0) 
(or equivalently between {q = 0,u = 1) and {q = l,u = 1)), as we see from Tab. [TJ 
is proportional to p = T/M. In this case, as it is remarked in the last paragraph, the 
fluctuations of the fleld q are getting stronger as the temperature increases, and this can help 
the system to jump from one vacuum to another and may have as a result the restoration 
of the Zm{2) symmetry (then the system passes from region (A) to region (B) in Fig. Q]). 

^We remind readers, that a 3 dimensional scalar field theory that includes powers of the scalar field up 
to four is supercnormalizable. If the powers of the scalar field are up to six the theory is renormalizable else 
the theory is nonrenormalizable, see for example Ref. jT,. 



5 The lattice model 

Our aim is to perform numerical computation and to see if the model can confirm the basic 
features of the phase diagram of Fig. ^ The only way to use this model for numerical 
computations is to discretise the action of Eq. (31) in the lattice. 

We will denote the lattice points by n = {ni,n2,n^) where ni, 712,713 are integers. If 
the lattice spacing will be denoted by a, then the corresponding physical points are Xn = 
{aril, 077.2, 0773). We remind readers that in this model we assume fixed lattice spacing which 
is set equal to the inverse value of the momentum cut-off A of the five dimensional gauge 
theory. 

The lattice action reads: 

Seff[q, u\= (3g -7r^^\ -{q{xn + ae^) - q{xn)f + p {u{xn + ae^) - u{xn)f 

+ H Veff{q{Xn),u{Xn)) (33) 

n 

where [3g = '^a/gl, e^ are the unit vectors, and we choose to measure all the quantities in 
the action in units of the lattice spacing. 

The potential Veff{q,u) is defined as Veff{q,u) = i^Veff{q,u). According to this defi- 
nition we have: 



where 



K//(g, u) = K7/(9) + Kffiu) + V:ff{q, u) (34) 

Ky/g) =4(rf-2)777'^~V457?^Esin^(vrrg)[l + E (^2 + ,;/^.)./. ] (35) 

VCS(. .) ^ -ie(. - 2),„-p^ 1 1 "ffyg:^' (37) 

Note that we have set t = Ta and m = Ma. 

The order parameters, according to which the symmetries Zt{2) and ^m(2) are violated 
or not, are defined by the average values of the Polyakov loops (see Eqs. (7) and (8)) which, 
for the lattice model we examine, are given by the equations: 

Pr{x) = cos{nu{xj), Py{x) = cos(7rg(x)) (38) 

The volume averages of the Polyakov loop are: 



Pr = ^T.Prix^n), Py = ^T. ^vM (39) 



and Vl = V/a is the lattice volume. 
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So with the lattice simulation we will measure the quantities: 

|P.| = (|P.|), \Py\ = {\Py\) 



(40) 



The critical values of t (or critical values of m) are specified as the values which maximize 
the susceptibilities: 



XiPr) = Vl {{Prf) - {{\Pr\)f , X{Py) = ^L ((P.)') " i{\Py\)) 



(41) 



6 Numerical results 




Figure 3: The phase diagram of the model for fixed l3g 
Vr = 163. 



0.2 or 0.1, for finite lattice volume 



In this section we will study numerically the phase diagram t—?Ti, for fixed fig. For this we 
have performed lattice simulation for several lattice volumes Vl = 10'^, 12^, 16^, 20^ and 24^. 
Near the peaks of the susceptibilities, for Vl = 24^, where the phase transition happens, we 
have used samples of 250K measurements which are separated by nine Metropolis iterations. 
The first 20K measurements were ignored for thermalization. For the other lattice volumes 
we have used samples with fewer measurements. 

In Fig. Elwe have plotted the phase diagram for (3g = 0.2 and 0.1 for finite lattice volume 
Vl = 16^. We see that it exhibits the main features of Fig. ^ namely it separates the t — m 
plane into the four distinct regions that correspond to the four phases (A),(B),(C) and (D) 
of the theory. 

The order parameter \Py\ as a function of m, for (3g = 0.2 and t = 0.5 is shown in Fig. |3] 
We observe that the phase transition of the system remains continuous even for large volumes 
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Figure 4: \Py\ and |Pt-| as a function of m for /5g = 0.2, t = 0.5 and Vl = 12^,24^. We 
see that as the lattice volume increases the phase transition, which corresponds to \Py\ , is 
getting more sharp. The phase transition for |Pt-| is not shown in this plot. 

Vl = 24^. Note that we have performed computations for several other characteristic values 
of t and jSg and the continuous behavior of the order parameter as a function of m is the 
same. 

The corresponding susceptibilities are shown in Fig. El The location of the points of the 
phase transition and the estimation of the relative errors have been done by computing the 
susceptibility several times in the range where the peak is expected. However we have not 
used a histogram method as it is not applicable to the model we examine. 

From the behavior of the peaks of the susceptibilities it seems that we have a second order 
phase transition. According to the theory of finite size scaling (see for example Ref. 0) we 
expect that the peak depends on the lattice volume as x{Py)max = cVl for large values of 
Vl, where h = j/Su (for the definition of the critical exponents 7 and u see Ref. jH]). This 
behavior was confirmed numerically, for the model we examine, and it is shown in Fig. IHl 
The errors in the figure were estimated by the Jackknife method. The critical exponent b 
was determined by a linear fit and it was found to be, for (3g = 0.2 , b = 0.70 ± 0.04 for 
t = 0.5, b = 0.68 ± 0.04 for t = 0.25 and b = 0.64 ± 0.04 for t = 0.05 using the three bigger 
volumes. For f3g = 0.1 and t = 0.5 we found that b = 0.66 ± 0.04. These numerical values 
for 6 (0 < 6 < 1) indicate that we do not have a first order phase transition and possibly the 
phase transition is of second order. Moreover, this is a strong indication that these phase 
transitions belong to the same universality class, as the numerical values of these critical 
exponents are very close and lie in the range of errors. Finally, we note, that these values 
for b could be equal to the corresponding critical exponent of the 3d ising model b = 0.657. 
This numerical value for the critical exponent b was obtained in Ref. [TTj . 
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Figure 5: x{Py) as a function of m for /5g = 0.2, t = 0.5 and Vl = 16^, 20^ and 24^. As the 
lattice volume increases rric is moving slowly towards the left and, for large lattice volume 
Vl, it seems to tend to an asymptotic value. 

We argue also that the qualitative features of the phase diagram, of Fig. 01 can not be 
just finite size effects. For this we have plotted in Fig. [7| the critical value rric for several 
values of t = 0.05, 0.25 and 0.5, for Pg = 0.2 and 0.1, as a function of lattice volume Vl- We 
see that there are small displacements for rric but the arrangement of the critical values does 
not change as the the lattice volume increases. This indicates that, in the infinite lattice 
volume limit, the qualitative features of the phase diagram are preserved. 

Note that we have not used the data points in Fig. [7|in order to to determine the critical 



exponent, as a fitting of the form rric 



rrir 



>+c7v^ 



l/3u 



will give unreliable results. The reason 



is that we have to determine three independent parameters whose values are very sensitive 
and we have not enough data points with a satisfactory accuracy. 

7 Conclusions 

We have computed the one loop effective potential for a 5D SU(2) gauge field theory at finite 
temperature and radius in the case of a background field with two constant components Ay 
and A^. 

However the effective potential, which is a perturbative result, can not explain straight- 
forwardly all the qualitative features of the phase diagram of Fig. [T] For this we constructed 
a phenomenological model by adding to the one loop effective potential kinetic terms (see 
Eq. (31)). We performed Monte Carlo simulations with the lattice version of this model (see 
Eq. (33)) and we found a phase diagram, for fixed lattice spacing and jSg, which exhibits all 
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Figure 6: Ln{x{Py)max) as a function of LniVL). The discrete points corresponds to Vl = 
10^, 12^, 16'^, 20^ and 24^. The continue hne corresponds to a best fit curve of the form 
Ln{xmax) = a + b ■ LniVL) for the three bigger volumes. 

the quahtative features of Fig. [T] 

Now the restoration of ^Af(2) for large temperatures (or the passing from region (A) 
to region (B) in Fig. Q) can be understood in the following simplistic way: even though 
the barrier that separates the vacua of the dimensionless field q increases linearly with the 
temperature T, the fluctuations of g are getting stronger (see section 4), and as it is confirmed 
by the lattice model, it succeeds in restoring the Zm{'^) symmetry. At the same time the 
fluctuations of u are getting more and more restricted so the field u is frozen to one of its 
vacuum states. 

Finally we remark that the numerical results indicate second order phase transitions, 
which is an interesting feature of the lattice model of Eq. (33). However, questions like the 
continuous limit (or nonperturbative renormalizability) are beyond the scope of this paper. 
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Figure 7: m^ for fixed f3g = 0.2, t = 0.5,0.25,0.05 and (3g = 0.1, t = 0.5 for five distinct 
values of lattice volume 10^, 12^, 16^, 20^ and 24^. The lines are there to guide the eye. 

9 Appendix: The constraint effective potential for two 
scalar fields 

Generalizing the definition of the constraint effective potential V^ff in Ref. nTH for the case 
of two order parameters we have: 



-VPRVeff{t4,t5) 



where 



and 



1 
1 /•« 



VA'^JiU - U)Kh - h)e 



-S[Al 



i^(^AX) = / dy I dPx TrPe*9^/o'^*("'"'2^)'^" 



UA 



2VR 

1 r^ 



I dr fd^x TrPe*3^r^5(x,r,?/)d?/ 



2Vf5Jo 
If we use the following representations for the delta functions 

r+oo 



and 



we obtain 



5{U -U)= I rfA4e'^^(**-*^) 

J —oo 

/+00 _ 

-oo 



(42) 

(43) 
(44) 

(45) 
(46) 

(47) 
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In order to compute the above path integral we use the saddle point method. We split the 
fields into classical and quantum parts 

^" = i?M + Ql (48) 

and 

A4,5 = &4,5 + g4,5 (49) 

We expand the exponent in Eq. (47) up to quadratic terms in quantum fields. The linear 
terms which are proportional to the equations of motion, according to the saddle point 
method are required to vanish (for details see Ref. ^U] )• 

^^4,5 T-r^ = (50) 



^ 6S 



. c. / =0 (51) 

i4,5-t4,5(5r,5)=0 (52) 

We assume that the background fields B^^ are constant and have different directions in the 
isospin space (note that we have also assumed that -B°2 3 = )• Of course there is not a 
gauge transformation that can put the two gauge field components in the same direction. 
However we can perform two rotations in the isospin space. With the first we can put B4 
towards the generator T3 and with the second we can put B^ on the plane defined by the 
generators T3 and Ti. Thus we can write 

^4 = BIT3 (53) 

and 

B, = Bin + BlT, (54) 

We will look for a saddle point solution assuming that it is possible to have 645 = 0. From 
Eqs. (50) and (51) we see that the classical field should satisfy the equations of motion for 
the gauge fields: 

- (^) ^ = (^^^- + 9^^''''K^K),.^ = (55) 

or 

Taking into account Eqs. (53) and (54) the only nonzero component of the field tensor is 

Fl = g,BlBl (57) 

From Eqs. (56) and (57), for /i = 5 and a = 1, we have 

{BlfBl = (58) 

and, for /x = 4 and a = 3, 1, we have 

(Bl^Bl = 0, a = 3 (59) 
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and 

BlBlBl = 0, a = 1 (60) 

These three equations are satisfied only if i?| = or 5^ = 0. However in general B\ ^ 
(If we assume that i?| = then we have only one gauge field component and the case is 
trivial). So we must have i^g = and thus the two gauge field components B^ and B^ have 
the same direction in the group space. 

Thus the saddle point we obtain has the form 



^4 = BiT^, B, = Bin, 64,5 = (61) 

Now if we compute the path integral of Eq. (47), around this saddle point, to one loop 
order, according to Ref. (TUI, we see that the result for the constraint effective potential is 
in agreement with that of Eq. (24) for the one loop effective potential. 
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